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ABSTRACT 

We analyze images of BIMA 12 C0 (J = l->0), VLA H I, and Spitzer 3.6 and 24 jxm emission 
toward the edge-on galaxy NGC 891 and derive the radial and vertical distributions of gas and the 
radial distributions of stellar mass and recent star formation. We describe our method of deriving 
radial profiles for edge-on galaxies, assuming circular motion, and verify basic relationships between 
star formation rate and gas and stellar content, and between the molecular-to-atomic ratio and hydro- 
static midplane pressure, that have been found in other galaxy samples. The Schmidt law index we 
find for the total gas (H2 + H I) is 0.85±0.55, but the Schmidt law provides a poor description of the 
SFR in comparison to a model that includes the influence of the stellar disk. Using our measurements 
of the thickness of the gas disk and the assumption of hydrostatic equilibrium, we estimate volume 
densities and pressures as a function of radius and height in order to test the importance of pressure 
in controlling the ph 2 / Phi ratio. The gas pressure in two dimensions P(r, z) using constant velocity 
dispersion does not seem to correlate with the ph 2 /phi ratio, but the pressure using varying velocity 
dispersion appears to correlate with the ratio. We test the importance of gravitational instability 
in determining the sites of massive star formation, and find that the Q parameter using a radially 
varying gas velocity dispersion is consistent with self-regulation (Q ~ 1) over a large part of the disk. 

Subject headings: galaxies: individual (NGC 891) — galaxies: ISM — galaxies: kinematics and dy- 
namics — stars: formation 



1. INTRODUCTION 

Edge-on galaxies have long been recognized as powerful 
probes of galaxy formation and evolution. The thickness 
of the s tellar disk is sensitive to the merging history of a 
galaxy (|Yoachim fc Dalcanto n 2006) , while the thickness 
of the interstellar medium ( ISM) may reveal the im print 
of gravitational instability (D alcanton et al.l l200T). The 
vertical structure of galaxy bulges can als o plac e con- 
straints on their formation: iCombes et al.l (|1990D have 
used N-body simulations to show that "boxy/peanut" 
bulges seen in many edge-on spiral galaxies may form 
from disk material via a bar instability. 

There are also practical advantages to observing edge- 
on systems. One can often achieve a higher signal to noise 
ratio (S/N) for detecting disk emission because more sig- 
nal is integrated along the line of sight. While this may 
be obvious for optically thin emission, even for H I and 
CO the optical depth may be moderated by large line-of- 
sight velocity gradients. Second, edge-on galaxies are the 
best objects to explore the vertical scale heights of disks, 
which are difficult to study in face-on galaxies. Determin- 
ing the mass surface density and the scale height allows 
one to derive average mass volume densities, essential for 
modeling galactic dynamics and gravitational collapse in 
the ISM. However, it is true that, obtaining radial dis- 
tributions of edge-on galaxies is not simple, since many 
radii contribute to each line of sight. Moreover, dust 



attenuation may adversely affect our ability to infer ra- 
dial profiles under the assumption of axisymmetry. In 
addition, a radial variation in the vertical distribution 
may not be apparent from an edge-on perspective. Use 
of kinematic information or deconvolution methods are 
therefore required. 

This is the first in a series of papers investigating how 
the gas layer thickness varies within edge-on galaxies and 
the resulting implications for star formation. In this pa- 
per we focus on NGC 891, a bright, nearby (D ~ 9.5 
Mpc fo r H = 75 km s" 1 M rxrMvan der Kruit fc Sear Id 
I1981H ). edge-on (i > 89° . lOosterloo et all 120071 ) spiral 
galaxy that has been studied extensively at various wave- 
lengths. We use H I images from the Very Large Ar- 
ray (VLA) and CO images from the Berkeley-Illinois- 
Maryland Association (BIMA) array to trace the atomic 
and molecular ISM, respectively. The existence of a ver- 
tically extended H I halo (a few kpc in extent) has been 
demonstrated by many authors (e.g.. | Sancisi & Allen 
119791: iSwaters et al.lll997t lOosterloo et al.ll2007| ). On the 
other hand, the vertical extent of the CO gas is still be- 
ing debated. While some authors have suggested that 
there is a thick compone n t to the CO gas lay er (e.g., 
iGarcia-Burillo et al.l Il992t iSofue fc Nakal I1993H others 
have argu ed that the CO gas layer is not extended but 
thin (e.g.. iScoville et alJll99l . We compare the verti- 
cal distributions of atomic and molecular gas by fitting 
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Gaussians to CO and H I intensity profiles taken in the 
vertical direction and obtain the gas disk thickness as 
a function of radius in order to estimate the gas volume 
density and pressure. These estimates are then compared 
with the midplane pressure derived from surface densi- 
ties alone, as is ty pically done for face-on galaxies (e.g., 
iLerov et all 1200^1 . 

A major motivation for obtaining direct estimates of 
gas volume densities is to provide better prescriptions 
for the star formation rate (SFR) in galaxies. Several 
authors have suggested that a power law relationship ex- 
ists bet ween the surface densities of gas (£ gas ) and SFR 
(S SF r) (|Schmidtlll959t lKennicutt]|1989l) The m ost com- 
monly used relation is due to iKennicuttl (|1998l) . who in- 
vestigated the star formation law (Schmidt law) in 61 
normal and 36 starburst galaxies, with £ gas and Esfr 
averaged over the disks, and determined a Schmidt law 
index of 1.4: Ssfr oc (S gas ) 1,4 . On the other hand, 
using spatially resolved data, iWong fc Blitd (|2002fi have 
sho wn that Ssfr is bet ter correlated with Eh 2 than Ehi, 
and ILerov et al.l ()2008f ) have reported a relationship be- 
tween star formation efficiency (SFE = £gFR/£ gas ) and 
stellar surface density, outside the "transition" radius 
where Eh 2 =Shi- These results indicate that a simple 
dependence of Ssfr on £ gas is an oversimplification, and 
that a more accurate description must take into account 
the role of the stellar disk's gravity in compressing H I 
gas to the high densities traced by CO. Consistent with 
this view, it has been noted that the radial distribution 
of i? m oi=SH 2 /^Hi correlates with the hydrostatic mid- 
plane pre ssure, which has contributions from both gas 
and stars (IWong fc Blitz|[2002T: iBlitaTfe Rosolowsky|2004t 
ILerov et al. 2008). We test these ideas using a Ssfr pro- 
file derived from 24 /im emission, a good SFR indicator 
(jCalzetti et al.1120071) . 

Moreover, star formation in disk galaxies appears to 
be suppressed beyond a threshold radius that is com- 
parable to the optical radius, even though the gas 
disk extends out much further. The physical origin 
of this star fo r mation threshold i s st ill widely debated 
(ISchavd l200l iBigiel et all l2010h . iKennicuttJ (fl989l ) 
and iMartin fc Kennicuttl (|2001| ) have suggested that the 
threshold depends on the gravitational instability of a 
galactic disk. They have investigated a relationship be- 
tween the star formation threshold and axisymmetric 
gravi tational instab ility based on the Toomre Q param- 
eter (jToomrel ll964) of the disk. The Q parameter for a 
thin rotating gas disk is defined by 



where n is the epicyclic frequency and a g is the gaseous 
velocity dispersion. Where the Q gas parameter is less 
than 1, instability is expected. We discuss whether the 
gravitational instability theory can explain the size of the 
CO disk or star formation disk in this paper. 

This paper is organized as follows. In £]2l we describe 
the CO and H I observations and data reduction and 
show maps of CO, H I, and Spitzer IR data. Section 
[3] shows vertically integrated position- velocity diagrams 
and the rotation curve. Sections 14.11 and 14.21 explain 
how we derive radial density distributions for the ISM 
from CO and H I data and for the stellar surface density 
(£*) and Ssfr from Spitzer 3.6 /im and 24 /jm data. In 



£15.11 we find the thicknesses of CO and H I and examine 
whether they have an extended thick component. In Sec- 
tion 15.21 and 15.31 we determine disk thicknesses and ver- 
tical velocity dispersions of gas and stars as a function of 
radius. Section 16.11 shows how the interstellar midplane 
pressure correlates with the molecular to atomic gas ra- 
tio. In £16.21 we examine relationships between Ssfr and 
HI, H2, and total gas, and compare a theoretical esti- 
mate of the star formation rate with the derived SFR. 
In £16.31 we show how the gravitational instability varies 
with the galactic radius in different circumstances: con- 
stant and varying velocity dispersions of gas and stars. 
In £16.41 we investigate how the interstellar pressure is 
related to the ratio Ph 2 /Phi in the vertical direction. 
Finally, we discuss and summarize our results in §£|7]and 
[HI respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

The 12 CO (J = l->0) observations of NGC 891 
were carried out using the 10-element Berkeley-Illinois- 
Maryland Association (BIMA) interferometer in 2002. A 
six pointing mosaic extending from near the center (a — 
2 h 22 m 33?3 and S = 42°20'52" at J2000) along the south- 
ern half of the galaxy (PA = 23°) was done in 3 configu- 
rations (B, C, D). The heliocentric systemic velocity V sys 
is 530 km s" 1 . 0136+478 (9.7° away, 3.8 Jy), 0359+509 
(19° away, 5.6 Jy), and Mars were observed as the phase, 
passband, and flux calibrator, respectively. The CO data 
were reduced using the MIRIAD package. The achieved 
channel maps (cube) from the MIRIAD task INVERT 
have angular resolution of 7" x 7" using natural weight- 
ing with 1" pixel size and 10 km s _1 velocity resolution. 
Figure [T] (top panel) shows the integrated intensity map 
rotated 67° (= 90°— PA) counter-clockwise and includ- 
ing channels from 220 km s _1 to 850 km s _1 . We use 
x and z to denote offsets parallel and perpendicular to 
the major axis of the galaxy (the galactic center, a = 
2 h 22 m 33?7 and 5 = 42°20'54", is placed at x = and 
2 = 0). The southern disk is placed in positive x. Since 
we adopt D = 9.5 Mpc, 1" corresponds to about 46 pc. 
In order to reduce noise in the intensity map, we used a 
masking method which blanks regions that fall below a 
3ct threshold in a smoothed (to 15" resolution) version 
of the cube. 

Although we lack a single-dish CO map of sufficient 
sensitivity to compare with the BIMA map, we have 
compared the integrated intensities along the major axis 
with the single-dish meas urements from the IRAJVQ 30- 
m telescope published by iGarcia-Burillo et "all (|1992f ) to 
check for flux being resolved out by the interferometer. 
We used the Dexter tool to extract data points from the 
published paper. The BIMA data were convolved to the 
IRAM beam size and sampl ed to match the positions ob - 
served by IRAM in Fig. 3 of lGarcia-Burillo et al.l (fl992l) . 
In fact, the interferometer flux (BIMA) agrees well with 
the single-dish flux over most of the southern disk, as 
shown in Figure [U The total BIMA flux in the figure is 
-0.9 times the IRAM flux. 

The H I data in B, C, and D configurations were ob- 
tained from t he NRAO Very Large Array (VLA) archive 
(|Rupenlll99lT ). We reduced the data using the AIPS (As- 
tronomical Image Processing System) and the MIRIAD 

1 Institut de Radioastronomie Millimctriquc. 
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X— offset (arcsec) 

Fig. 1 — Top: CO integrated intensity map of NGC 891. Contours are 20, 34, 59, 101, 174 and 298 K km s _1 in logarithmic scale. The 
lowest contour level is ~ 2a. Note that the positive x-offset values are to the south and the negative values are to the north. The dashed 
line shows the 50% sensitivity contour. The synthesized beam (7" X 7" ) is shown in the lower right corner. Middle: Zoom-out version 
of the top panel shown in the same scale with the H I map. The box indicates the region shown in the top panel. Bottom: H I integrated 
intensity map. Contours are 1.30 (~ 2a), 1.83, 2.58, 3.64, 5.14 and 7.25 X 10 3 K km s _1 . The absorption by SN 1986J is visible near 
x = 60". The synthesized beam (11"56 X 8'.'78) is shown in the lower right corner. 
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Fig. 2. — Comparison of integrated intensity along the major axis 
be tween our data (BIMA) and the single-dish data (IRAM) given 
bv lGarcia-Burillo et all 1199211 . 



packages. Using AIPS we calibrated the data and sub- 
tracted continuum emission deduced from 3 edge chan- 
nels without line emission. After the calibration and 
the continuum subtraction, the data were self-calibrated 
with the MIRIAD task SELFCAL and mapped with the 
MIRIAD task INVERT. In order to obtain a reasonable 
angular resolution without significant loss of sensitivity, 
we used a robust weighting with a robustness factor of 
0.4 and with a 2" cell size so the achieved beam size is 
ll'/56 x 8'.' 78. The resulting cube consists of 31 chan- 
nels from 220 to 840 km s _1 with a velocity resolution of 
20 km s -1 . The velocity integrated intensity map of the 
H I is shown in Figure [T] (bottom panel). As with CO, 
this H I map is the result of masking using a smoothed 
version of the cube, although the smoothing is at 20" res- 
olution in this case. The absorptio n at around x=60" is 
because of the supernova SN1986J (|Rupen et al.l[l987h . 

We have obtained near/mid-IR images at 3.6 /zm 
(IRAC) and 24 fj,m (MIPS) from the Spitzer archive: 
Program ID 3 (PI: G. Fazio) for IRAC and Program ID 
20528 (PI: C. Martin) for MIPS. Maps of Spitzer data 
are shown in Figure [3] When we derive the radial distri- 
butions of stellar mass (3.6 /im) and recent star forma- 
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Fig. 3. — 3.6 fim (top) and 24 (im (bottom) emission from NGC 891 as imaged by Spitzer. Contour levels are 0.30 (~ 5ct), 0.83, 2.29, 
6.31, 17.38 and 47.86 MJy sr _1 in logarithmic scale. The box in the upper panel indicates the region shown in the top panel (CO) of Fig. 

m 



tion (24 /im), foreground bright stars are first blanked 
and filled by interpolations of adjacent data points us- 
ing GIPSY (Groningen Image Processing System) tasks 
BLOT and PATCH, respectively 

3. KINEMATICS 

3.1. Position-Velocity Diagrams 

In order to study the kinematics of the galaxy, we 
have integrated the data cubes in the z dimension to 
derive vertically integrated position-velocity (PV) dia- 
grams. The PV diagrams in CO (integrated over ± 
10" from the plane) and H I (integrated over ± 20") 
are shown in Figure 01 The PV diagram of H I in the 
figure displays an asymmetry that indicates a more ex- 
tended disk on the southern side . As previously noted by 
IGarcia-Burillo et all (|l99l and lSofue fc NakaJ (fl99l . a 
central nuclear feature (fast-rotating disk) followed by a 
gap and a second peak is apparent in the PV diagram of 
CO. IAthanassoula fc Bureau! f|1999f ) have suggested that 
the presence of such a gap is evidence of a bar, by show- 
ing numerically simulated PV diagrams of edge-on barred 
galaxies seen at different viewing angles. The gap de- 
velops as gas is depleted from the outer bar region by 
gravitational torques from the bar. In addition, they 
have shown that a side-on bar generates higher velocities 
in the central emission feature than in the outer parts, 
while an end-on bar does not produce such a feature. 
Based on their results, we suggest that NGC 891 has a 
bar seen side-on rather than end-on. 

3.2. Rotation Curve 

We derived the rotation curve shown in Figure [5] using 
PV diagrams of CO and H I along the major axis. Note 
that unlike the vertically integrated PV diagrams shown 
in Fig. |4j these PV diagrams are slices along the mid- 
plane, in order to prev ent contamination b y more slowly 
rotating halo gas (e.g. JSwaters et aLlll997l) . The curve is 
the result of combining CO (red circles and dashed line) 
for the inner region and H I (blue crosses and solid line) 
for the outer region because the CO emission is strong 



near the center but weak in the outer part, while the H I 
emission does not have a prominent central component 
as seen in the CO data (but see discussion in Section 
|4~T|) . Only the southern part (V r > V sys ) of the PV di- 
agrams is used for obtaining the rotation curve, due to 
the coverage of our observations. The rot ation curve i s 
obtained by the envelope tracing method (|Sofud [i996'). 
which is based on a terminal velocity (Vt er ) corrected by 
the observational velocity resolution a b s (10 km s _1 for 
CO and 20 km s _1 for H I) and velocity dispersion of the 
gas (assumed to be a g = 8 km s _1 ): 

Vrot = Vter ~ + < (2) 

The highest-velocity 3<r contour in the PV diagrams is 
selected as the terminal velocity Vter- The rotation ve- 
locity appears to rise rapidly to a maximum velocity of 
~255 km s _1 at 10" (suggesting solid-body rotation of 
a nuclear disk), then decreases to a minimum velocity 
of ^180 km s -1 , followed by an increase again to a sec- 
ond peak. Beyond the second peak, the rotation curve 
flattens at about 2 30 km s -1 . The size of the correction 
term {J o- 2 ohs + aj) is used for representative error bars 

for CO (upper left corner in Fig. [5} and H I (upper right). 
Our adopted rotation curve, which is a combination of 
the CO and H I rotation curves, is used to determine 
the instability parameter Q (Section l6.3|) . although our 
results are not sensitive to the shape of this curve. In 
light of the discussion in Section 13.11 the rise and fall of 
the curve in the central region of the galaxy is likely due 
to bar-induced gas streaming motions rather than a real 
change in the mean circular velocity. 

4. RADIAL DISTRIBUTION 

As discussed in Section [TJ radial distributions are not 
trivial to infer for an edge-on galaxy. Two basic ap- 
proaches can be followed. If line-of-sight velocities are 
measured, emission at a particular velocity can be as- 
signed to a particular position along the line of sight 
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X— offset (arcsec) 

Fig. 4. — Top: CO position-velocity diagram integrated over the minor axis (± 10"). Contour levels are 0.030, 0.046, 0.070, 0.107, 0.163 
and 0.249 Jy arcsec -1 in logarithmic scale. The lowest level is ~ 3cr. The systemic velocity V sy a = 530 km s" 1 is shown by a horizontal 
dotted line. The vertical dashed lines show the 50% sensitivity. The box represents the excluded region explained in Section 14.11 Bottom: 
H I position-velocity diagram integrated over the minor axis (± 20"). Contours are 0.180 (~ 3c), 0.297, 0.489, 0.805, 1.326 and 2.185 mjy 
arcsec -1 in logarithmic scale. The absorption at around 60" is due to SN 1986J. CO contours (red) are overlaid on H I contours (blu e). 
The vertical dashed lines show the 50% sensitivity of the CO mosaic. The box represents the excluded region explained in Section UT] 
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Fig. 5. — Observational rotation curve from CO and H I. The 
red circles and blue crosses show the CO and H I rotation velocity, 
respectively. The representative two sided error bars are shown in 
upper left (CO) and upper right (H I) corners. 

(given a model of the galaxy's rotation), and thus a par- 
ticular galactocentric radius. If no velocity information 
is available, an inversion technique can be used to derive 
the radial profile from the pr ojected brightne ss distribu- 
tion, assuming axisymmetry ( Warmcls 1988). We apply 
the first method to the CO and H I data and the second 
method to the 3.6//m and 24/im data. In Appendix [A] 



we perform comparisons of the two methods on actual 
and simulated data. Due to limitations inherent in both 
methods, we estimate that individual points in our de- 
rived radial distributions are likely to be accurate to only 
a factor of ~2. 

4.1. Molecular and Atomic Gas 

We derived the radial gas distributions using the 
position-velocity (PV) diagrams. We hereafter refer to 
the method using the PV diagram as the PVD method. 
The PVD method makes use of more data than taking 
a strip integral near the te rminal velocities (e.g. iRupenl 
119911 : ISofue fc Nakailll993f ) which discards data at inter- 
mediate velocities which still contribute to the flux of 
the galaxy. For comparison with the atomic gas (H I) 
profile, the CO data have been convolved to the H I 
beam size (ll'/56 x 8'/ 78). The PV diagrams have been 
produced by integrating the data cube in the vertical 
(minor axis) direction without masking (±10" and ±20" 
from the plane for CO and H I, respectively) and the 
PVD method assumes circular rotation and a flat rota- 
tion curve which is a reasonable approximation based 
on Figure [5] Each pixel in the PV diagram can be as- 
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Fig. 6. — (a) Radial profiles of H2 (red open squares), H I (blue open circles) and total gas (solid circles) surface density, at the resolution 
of the H I data. Horizontal and vertical error bars for the H2 and H I profiles are explained in the text. The vertical error bar in the 
lower right corner represents the adopted uncertainty for the total gas, based on the difference between PVD and RADPROF profiles (see 
Appendix0- The solid lines (red for CO and blue for H I) represent the 3-sigma detection threshold, (b) Stellar surface density obtained 
from the task RADPROF (blue solid line) and from the exponential disk model (red dotted line) as a function of radius based on the 3.6 
(im emission. The bulge part is excluded in the model fitting. The error bar in the lower left represents the uncertainty, (c) SFR surface 
density as a function of radius based on the 24 fim emission. The representative error bar in the lower left corner indicates a factor of 2 
change obtained from the biggest difference between RADPROF and ELLINT profiles for several 24 (im maps of face-on galaxies. 



sociated with a galactocentric radius using the observed 
radial velocity (V r ) and the assumed circular speed (V c = 
250 km s^ 1 ) at each position x: 

r = Vc ( v -v ) with 1 Vr ~ Vs ^ < Fc ' (3) 

\ *r * sys / 

where the mean value of xj (V r — V sys ) within a pixel with 
velocity width AV r (10 km s" 1 for CO and 20 km s" 1 
for H I) can be estimated as 

/ X \ = ]x\_ , / l^-^sy S |+AK/2 \ 

\ V r - V sys I AV r \\V r - V sys \ ~ AV r /2j ' 1 ' 

The flux of each pixel in the PV diagram is corrected to 
the corresponding surface brightness in a face-on galaxy 
after considering the depth along the line of sight and 
the fact that both near and far sides contribute to one 
pixel. In this procedure, a region near the center of the 
PV diagram (bounded by \x\ < 80" and \ V r - V sys \ < 80 
km s _1 ) shown in Figure 0] as a box is excluded due to 
blending of emission from many radii. 

We have derived H2 surface mass density from the in- 
ferred face-on surface brightness, usin g a conversion fac- 
tor d e rived from Galacti c observations ([Strong fc Mattoxl 
[19961 IDame etalll200l : 

N(U 2 ) [cm,- 2 ] = 2 x 10 20 Ico [Kkm s -1 ]. (5) 

To obtain H I surface density, we use the optically thin 
approximation : 

iV(HI)[cm" 2 ] = 1.82 x 10 18 Im [Kkm s -1 ]. (6) 

The radial distributions of H2, H I and total gas 
(H2+H I) are shown in Figure [Bfa). Each point rep- 
resents an averaged value of data in a 12" radial bin. 
Horizontal error bars in the figure show maximum and 
minimum values of radius as an uncertainty derived by 
varying x and V r in Equation §3§ by the angular (Ax) 
and the velocity (AV r ) resolutions. Vertical error bars of 
H2 and H I profiles reflect dispersions of the fluxes within 



each annulus (standard deviation). The solid lines (red 
for CO and blue for H I) represent 3ct detection thresh- 
olds obtained using the la noise of the PV diagrams as 
input into the PVD method; these become larger near 
the center because smaller line-of-sight depths associated 
with smaller x- values in a given velocity interval are used 
to normalize the fluxes. Therefore, a few data points of 
the H I profile near the center are considered unreliable 
because they are below the detection threshold. How- 
ever, inclusion of these data points have little effect on 
the total gas profile. 

Examination of the CO profile reveals an obvious 
concentration of molecular ga s toward the center fol- 
lowed by a gap , consi stent with lScoville et all (|1993l ) and 
iSofue fc Nakail (|1993l ). and a second peak at a radius of 
around 4.6 kpc (100"). Beyond the peak, it decreases 
slowly toward outer radii. CO extends to ~ 280" in 
radius; the gain falls to 50% at th is radius so the pro- 
file is truncated beyond this radius. iGarcia-Burillo et "all 
(1992) reported that the "broad ring-like structure" at a 
distance of 3.6 to 4.0 kpc from the galact i c cent er may 
be caused by spiral arms. ISofue fc Nakail (fl99l inter- 
preted the secondary peak as a "3.5-kpc molecular ring" 
(with D = 8.9 Mpc) and noted several intensity peaks 
beyond the molecul ar ring, possib l y re sulting from spiral 
arms. In addition, IScoville et al.l (|1993t ) have suggested 
the peak at 4 kpc may be the result of "a spiral arm 
tangential to the line of sight" . 

The atomic gas profile also shows a gap corresponding 
to the gap in the CO profile. There appears to be some 
H I near the center, although the errors are larger for the 
reaso ns discussed above. While sever al earlier studies 
(e.g.. lRupenll!991t ISofue fc N akai 1993) suggested a lack 
of atomic gas near the center, H I associated with the 
central CO disk can be seen in the Westerbork Synthe- 
sis Radio Telescope (WSRT ) PV diagrams pres ented by 
IS waters et all (|1997l ) and lOosterloo et all (|2007l) . Unlike 
the molecular profile, the atomic distribution outside of 
the gap stays nearly flat until it decreases in the outer 
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Fig. 7. — Top: Exponential disk model fvan ■ 
2.88 and 18.20 in unit of MJy sr" 1 . Middle: Residual map obtained by subtracting the model from 3.6 /im data. Contour levels are same 
to the top panel. Bottom: Residual map with contours of negative values: -0.18, -0.46, -1.15, -2.88. 



disk. 

The total gas surface mass density is estimated by com- 
bining the H2 and H I profiles and including associated 
helium: E gas = 1.36(Sh 2 +^hi)- Since the CO profile has 
been obtained from the southern disk (the northern disk 
is unavailable in our data), we used only the southern 
disk of H I when deriving the radial profile to compare 
with the CO profile. Consequently, the total gas profile 
is obtained only for the southern half of the galaxy. As 
shown in Figure a), the total gas distribution closely 
tracks the molecular gas distribution, especially within 
r=200". We use a factor of 2 as the representative error 
for the total gas profile, based on our tests described in 
Appendix |A1 This does not include possible errors in our 
conversion factors (Eq. [5] and [6]). 

4.2. Infrared: 3.6 and 24 \xm 

In order to obtain radial density distributions from the 
Spitzer data, which lack kinematic information, we used 
the GIPSY task RADPROF. RADPROF takes an inte- 
grated intensity strip for both the left and right side of a 
galaxy (if available) and compute s a radial d is tribution 
using the Lucy iterative scheme (|Lucvl fl974l ; iWarmelsl 
1988). There are four initial functions that can be used 
in RADPROF: linear regression, exponential decreasing, 
Gaussian distribution, or flat distribution. In our anal- 
ysis, an exponential decreasing function for 3.6 /im and 
the Gaussian distribution for 24 /im were used as initial 
guesses before iterating. Our results are not sensitive to 
this choice. The task assumes an axisymmetric disk. Be- 
fore employing RADPROF, the 3.6 (FWHM = l'/6) and 
24 /im (FWHM = 5'/9) data have been convolved to the 
H I beam. 

We have also used an isothermal and s elf-gravitating 
disk model (jvan der Kruit fc Searldfl981af) to derive the 
radial distribution of 3.6 /im emission by fitting an expo- 
nential disk model. The central part, within ±50" from 
the center, is excluded in the fitting due to the presence 
of a stellar bulge. The fitting function for an edge-on 



galaxy, 

Im(x, z) = 0(0, 0) (|) K X (£) sech 2 (£) , (7) 

is obtained by integrating the model, 

L(r, z) = L e- r/l sech 2 (£\ (8) 

along the line of sight, where Lq is the space luminosity 
density at the center, I is the scale length in the radial 
direction, z* is the scale height independent of radius r, 
/i(0, 0) = 2hL , and Ki is the modified Bessel function 
of order 1. The scale length (I) and height (z*) obtained 
from fitting the 3.6 /im data are about 80" (~ 3700 pc) 
and 8" (~ 370 pc), respectively. Fitting the projected 
intensities (Eq. [7]) of an exponential disk model (Eq. 
[8|) to the 3.6 um image, we made model and residual 
maps as shown in Figure [7) The residual map shows an 
asymmetric thick component skewed towards negative x, 
which might represent the near side of the bar. In the 
inner disk, the fit is affected by the bulge light, creating 
positive residuals at high z and negative residuals at low 
z. In the outer disk, the assumption of constant stellar 
scale height means that the model disk is thinner than 
the actual disk (see Section I5T21 . 

Figure [l^b) shows the derived stellar surface mass den- 
sity profiles. The exponential disk model is shown as a 
dotted line and the RADPROF solution as a solid line. 
Their profiles agree well in the disk. The stellar ra- 
dial profile from the exponential disk model will be used 
throughout this paper. The representative error bar in 
the lower left corner is obtained from the biggest differ- 
ence between the two profiles in the disk, about a factor 
of 1.5. We convert from the units of the Spitzer image 
(MJy sr _1 ) to stellar mass density (£*) using a conver- 
sio n factor for 3.6 /im intensity (I3.6) empirically derived 
bv lLerov etall (|2008D : 

[M pc" 2 ] = 280 (cos i) I 3 . 6 [MJy sr" 1 ], (9) 



where the inclination of the galaxy i = 0° since we have 
obtained the radial density profiles using the exponen- 
tial model and RADPROF, which converts an edge-on to 
face-on. In doing so . we assume 3.6 /im is dominated by 
old stars. Use of the lLerov et al.l (|2008D conversion factor 
is supported by its consistency (to within 30%) with an 
independe nt, J — K color base d mass-to-light ratio em- 
ployed by Ide Blok et al.l (|2008f ) . We have compared the 
median 3.6 /im intensity with the median if -ba n d (2.2 
/im) intensity in 10" bins following ILerov et al.l (|2008[) 
and obtained a I^.q/Ik ratio of 0.68. ILerov et al.l (2008) 
used a value of Is.s/Ik = 0.55 to convert from 3.6 /im 
intensity to stellar mass density £*. Even though there 
is a discrepancy between two ratios (a factor of ~1.2), 
our data show a linear relationship between 3.6 fim and 
-ftT-band intensities. Therefore, the errors introduced by 
ignoring dust attenuation in the near-infrared appear to 
be within the uncertainties (a factor of 1.5) assumed for 
t he 3.6 um radia l profi le . 

ICalzetti et al.l (|2007f ) have examined correlations be- 
tween mid-IR emission (8 and 24 /im) and star forma- 
tion rate (SFR) and concluded that mid-IR emission, 
especially 24 /im , is a good tracer of SFR. For the con- 
version from 24 /im luminosity surface density S^Aam 
(ergs -1 kpc -2 ) to SFR s urface density Ssfr, we adopt 
the calibration derived by ICalzetti et al.l (|2007l) : 

/ q \ 0.8104 

M 2 = 1-56 x 10-- ( %g \ , 

M.0 yr 1 kpc 1 \ erg s 1 kpc z ) 

(10) 

where 

and I v is the 24 /im surface brightness derived from the 
RADPROF solution. The derived S S fr (Eq. HUJ) as a 
function of radius is shown in Figure [SJc) . The profile 
of Esfr falls sharply from the center to around 40" and 
then rises slightly before decreasing more gradually in 
the outer disk. This tendency seems similar to the CO 
distribution rather than H I, implying that SFR is more 
strongly correlated with the molecular gas surface den- 
sity. We also derived another SF R radial profile using 
the 1.4 GHz radio continuum map (|Dahlem et al.|[l995f) 
to check whether the 24 /im SFR profile is consistent 
with it. Since Ha and UV data show severe extinction in 
edge-on galaxies, the radio continuum is perhaps the best 
alternative for comparing with the SFR inferred from the 
24 /im data. For obtaining the SFR from the radio con- 
ti nuum, we used the m ethod described in the literature 
bv lMurgia et~aTl (pOOl and the RADPROF method. We 
confirmed that the ratio of 24 /im to radio continuum for 
total SFR is about 1.1. 

Since the distribution of recent star formation in a 
galaxy may deviate strongly from axisymmetry, which 
is a basic assumption of the RADPROF method, we 
tested the ability of RADPROF to recover the azimuthal 
averaged radial profiles of several 24 /im maps of face- 
on galaxies: NGC 628, NGC 2403, NGC 3184, NGC 
4321, NGC 4736, NGC 5055, NGC 5194, NGC 6946 and 
NGC 7331. These galaxies have been observed as part 
of the Spitzer Infrared Nearby Galaxies Survey (SINGS) 
(|Kennicutt et al.l 12003) . We compared the radial pro- 
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Fig. 8. — Single Gaussian fitting to the CO and H I integrated 
maps. Three panels show the integrated intensity, vertical offset, 
and deconvolved FWHM thickness from top to bottom. The ver- 
tical dotted lines in the first panel show the 50% sensitivity limit 
for CO. The horizontal dashed lines (for H I) and dotted line (for 
CO) in the bottom panel represent the weighted mean values. 

file obtained from the GIPSY task ELLINT integrating 
a map in elliptical rings with the profile derived from 
the RADPROF solution using a strip integral. The two 
different methods yield quite consistent results, suggest- 
ing that asymmetries due to disk substructure do not 
prevent determination of the radial profile. The stan- 
dard deviation of the differences for a particular galaxy 
ranges from 0.1 dex (NGC 3184) to 0.3 dex (NGC 2430). 
We use the largest of these values (0.3 dex, or a factor of 
2) to represent the error in the 24 /im profile. 

5. VERTICAL DISTRIBUTION 

We now proceed to measure disk thicknesses and verti- 
cal velocity dispersions as functions of radius. The mea- 
sured thicknesses of the CO and H I disks enable us to 
infer volume densities directly, rather than from the as- 
sumption of hydrostatic equilibrium. The radial distri- 
bution of velocity dispersions for gas and stars is more 
difficult to determine observationally but can be calcu- 
lated for a couple of different model assumptions. We 
describe the models in this subsection and then use in 
Section [S] In addition, models with constant velocity 
dispersions of gas and stars will be included for compar- 
ison purposes. 

5.1. Disk Thickness in Integrated Intensity 

To establish the overall characteristics of the CO and 
H I layers, we have fitted single or double Gaussian pro- 
files along the vertical direction to the integrated inten- 
sity maps of CO and H I using the MIRIAD task GAU- 
FIT. For the CO and H I data, we used the masked 
maps shown in Figure [T]to obtain higher S/N in order 
to yield more reliable fits. Figure [8] presents the results 
of single Gaussian fitting to the CO and H I maps. The 
integrated flux, vertical offset (centroid of the disk), and 
deconvolved FWHM are plotted from top to bottom, re- 
spectively. The CO flux of the left side is falling sharply 
towards the outer disk due to lack of sensitivity in our 
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Fig. 9. — Left: Gaussian widths of H I (open circles) and H2 (filled circles) as functions of radius normalized by the optical radius 
''25 ~ 400". The dotted (H I), solid (H2) and dashed (total gas) lines show linear approximations obtained by least-squares fitting. The 
fitting functions are shown on the top. The shaded regions around the fitted lines represent uncertainties of their fitting. Right: Stellar 
scale height as a function of radius normalized by the optical radius. The red and blue solid lines show the scale heights of west and east 
sides, respectively. Average of the west and east disks is shown as filled circles. The dashed line indicates the linear approximation by 
least-squares fitting to the averaged data points and the shaded region around the best-fit represents uncertainty of the fitting. The fitting 
function is shown on the bottom. 

mosaic (50% sensitivity is shown as vertical dotted lines 
in the figure.) From the results of the single Gaussian 
fitting, we have obtained a weighted mean value for the 
FWHM thickness of - 7" (~ 320 pc) for CO and - 19" 
(~ 870 pc) for H I. Across the galaxy, the scatter around 
the mean value is about 1.3" (CO) and about 4.5" (H I). 
In the case of H I, the fit improves significantly if a sec- 
ond Gaussian is fit to the profile; the resulting FWHM 
thicknesses are around 14" (<~ 650 pc) for the thinner 
component and around 44" (~ 2 kpc) for the thicker com- 
ponent. The thinner component of the double Gaussian 
fit can be basically identified with the single Gaussian 
fit. The thicker component may be due in part to line of 



sight p rojection of a flaring disk, although iSwaters et al.l 



(1997) argue that an additional, slowly-rotating halo 
component is also needed. The vertical offs ets of the H I 
in the outer disk show evidence of a warp (jRupenl 119911: 
ISwaters et aTlll997l;IOosterloo et aLl l2007). We could not 
find clear evidence of an extended th ick CO disk as sug- 
;ested bv lGarcia-Burillo et al.l (|1992fl and lSofue fc Nakal 
1991 . since a double Gaussian profile is not needed to 
fit the CO image, implying a one-component thin disk 
layer. iScoville et al.l (|1993l ) also failed to find evidence 



for a vertically extended CO distribution. 



5.2. Radial Variation in Disk Thickness 

In order to determine the disk thicknesses of H2 and 
H I as functions of radius, we have derived the Gaus- 
sian width (0.42 times the FWHM) by fitting a Gaussian 
function to the z profile of a velocity-integrated intensity 
map. We have only integrated over the terminal veloc- 
ities on the redshifted side (740-760 km s _1 ) in order 
to provide a cleaner diagnostic of how the gas thickness 
varies with radius (rather than with x). In addition, 
fitting to the terminal velocities exclud es any slowly- 
rotating halo component (see Fig. 4 of Swat ers et al.l 



1997). However, a double-Gaussian fit to the terminal 
velocity maps proves unreliable due to limited signal-to- 
noise, so only a single Gaussian is fit, even for H I. Thus, 
we are effectively considering only the thin disk H I com- 
ponent in this section. More sensitive high-resolution 
H I data would enable a two-component fit of the verti- 
cal profile over this limited velocity range. The obtained 
Gaussian widths as functions of radius, normalized by 
the optical radius r 2 5 (« 400"), are shown in Figure|9ja), 
which indicates that both Gaussian widths increase mod- 
erately with radius. Each circle represents an averaged 
value of data in a 12" radial bin and the dotted (H I), 
solid (H2), and dashed (total gas) lines show linear ap- 
proximations obtained by least-squares fitting. The fit- 
ting functions are shown on the top of the figure. In order 
to obtain the scale height of the total gas, we used the 
combined map of CO and H I. Assuming that gas cycles 
between the atomic and molecular phases, the combina- 
tion might be considered as a single dynamic gas compo- 
nent. The minimum scale height of the gas is ~2.5" and 
the scale height at 0.5r25 is ~4.9". The shaded regions 
around the best-fit lines represent uncertainties in the 
fitting. Note that the Gaussian width is approximately 
2~ 1//2 times the scale height of a sech 2 fit. 

To estimate the variation in the stellar scale height 
with radius, we have fitted a sech 2 (z/z») function 
(|Spitzerlll942|) to the vertical (z) distribution of the 3.6 
(im data at each radius after obtaining radial profiles at 
different z. In order to obtain the radial profiles, we ran 
RADPROF for many different values of z from —33.6" to 
33.6" in steps of 2.4". The obtained scale height profiles 
for the west side (red line) , east side (blue line) and their 
average (filled circles) are shown in Figure HJb). Since 
the west and east profiles show differences that may be 
related to the bar structure, we use the average profile. 
After excluding the central region (r < 90") near the 
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TABLE 1 

Models of Vertical Velocity Dispersions 





Gas 

Value 


Model 


Stars 

Value 


Model 


Poisson 


<r s (r) given by Eq. 1121 


PG 


cr*(r) given by Eq. [12] 


PS 


Constant z 






<r„(r) = x/ttGS.z,, e» = 8" 


ZS 


Constant a 


u g = 8 km s 


CG 


cr* = 21 km s — 1 


cs 



References. — PG a nd P S: I Naravan fc Jo e (2002b); ZS: Bottcma (1993); CG: 
IShostak fc van der Kruit (1984); CS: Rafikovl {200% 

Note. — PG: Poisson Gas; PS: Poisson Stars; ZS: Constant z Stars; CG: Constant 
Gas; CS: Constant Stars 
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Fig. 10. — Vertical velocity dispersions obtained from Eq. I12l as 
a function of radius, normalized by the optical radius, for the gas 
(blue solid line) and stars (thick red solid line). The red dotted line 
shows stellar velocity dispersion using cr„ = y/nGH*z* IBottemal 
1993) with the constant scale height (8"). Each horizontal dashed 
line indicates a constant velocity dispersion for the gas (8 km s _1 ) 
and stars (21 km s _1 ). 

bulge, least-squares fitting to the averaged data points 
is used to obtain the linear approximation shown in the 
figure as the dashed line with the shaded region repre- 
senting the uncertainty of the best fit. The minimum 
scale height is 3.8"and it increases to ~15". 

5.3. Radial Variation in Vertical Velocity Dispersion 

A number of observations hav e shown that 
the veloci ty dispersions of gas (jPetric fc Rupenl 
120071: ITamburro et all 120091) and stars 

(jvan der Kruit fc Freeman! 1 19841 ) decline with ra- 
dius. We have derived vertical velocity dispersions as 
functions of radi us for the gas and stars using Equation 
(fl2j) provided bv lNaravan fc Jod (|2002bD : 



47rG?po ;to tPoi 
-(d 2 pi/dz 2 ) z=0 



(12) 



where the dark matter halo is ignored and boundary 
conditions at the midplane (pi = poi and dpi/dz = 0) 
are used. This equation implicitly assumes hydrostatic 
equilibrium between gravity and turbulent pressure. The 
subscript i can be replaced by either g or * for the total 
gas and stars, respectively. The total midplane density 



(Po.tot) includes p 0g [= H g / (h g V2n)] and p * (= £*/2z*), 
where h g and z* are the scale height of the gas and stars, 
respectively (see Appendix [Bj. As mentioned in Sec- 
tion l5.2l the gaseous and stellar volume densities follow a 
Gaussian distribution and a sech 2 function, respectively. 
The derived vertical velocity dispersions from Eq. [12] are 
shown in Figure H0l We refer to these models as PG and 
PS for gas and stars, respectively, since they are based 
on the Poisson equation. The gas velocity dispersion 
(blue solid line) falls off strongly with radius, contrary to 
the usual assumption o f constant a g . A recent study by 
ITamburro et al.l (|2009f ) also found that the H I velocity 
dispersion decreases as a function of radius. Our cen- 
tral velocity dispersions are comparable (within a factor 
of 2 or less) to values found in some of their galaxies. 
The thick red solid line represents the stellar velocity 
dispersion (<r» ) when a varying stellar scale height is ap- 
plied with Equation (|12l) . For comparison, we also show 
the stellar velocity dispersion assuming a constant stellar 
scale height (8" from Section |4~2|) in a purely stellar disk, 
so that er* = \J-kGY>*z* (|Bottemalll993[ ). This model is 
listed as ZS in Table HJ 

Many studies have assumed that the velocity disper- 
sions are constant with radius. So, we also plot the 
constant values of the velocity dispersions for compar- 
ison. Each horizontal dashed line represents a constant 
value of velocity dispersion for the gas (model CG) and 
stars (model CS). A constant value of <j q ~ 6-8 km s _1 
is widely adopted base d on observations of H I and 
CO in face-on galaxies (fShostak fc van der Kruitl fl984l : 
iCombes fc Becauaertj|1997l) . and for this model we adopt 
a value of o~ g = 8 km s _1 (providing only modest al- 
lowance for additional vertical support from cosmic rays 
or magnetic fields). For a constant value of stellar ve- 
locity dispersion, we adopt a value of cr* = 21 km s _1 
from th e upp er range (er*, r = 35 km s _1 ) quoted by 
I Rafikovl {2001) corr ected for the y elocity anisotropy us- 
ing tT*/<7* ir = 0.6 (Bottcmal ll993D . Models of velocity 
dispersions listed in Tabic Q] will be used in Section [6] 

Note that all models that we consider are based on ob- 
servations of face-on galaxies or on physical assumptions 
since we do not have direct measurements of the veloc- 
ity dispersions. A more direct approach t o deriving the 
gaseous velocity dispersion is discussed by I Oiling] (1996) 
and will be investigated in future work. 

6. THE H 2 /HI RATIO AND STAR FORMATION 

6.1. Hydrostatic Midplane Pressure 

.Previous studies (e.g..lElmegreenlll993l: iWong fc B litz 
I2002t iBlitz fc Rosolowsky||2004D have suggested that the 
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Fig. 11. — (a) Midplane pressure as a function of r adiu s normalized by the optical radius. The red solid and dashed lines sho w the 
hydrostatic midplane pressure obtained from Equation (1 1 3 P and the turbulent gas pressure at the midplane derived by Equation i21l using 
the constant gas velocity dispersion (Model CG), respectively. The black solid and dashed lines represent the hydrostatic pressure and the 
turbulent gas pressure when varyi ng v elocity dispersion (Model PG) is used. The blue circles show the hydrostatic pres sur e including the 
gas disk obtained from Equation (IB6I ) . (b) Ratio of Sh 2 /Shi as a function of the hydrostatic midplane pressure (Eq. 1136 with constant 
a g and 2*. The power law slope (dotted line) is 0.89. The vertical error bars indicate the uncertainties in the ratio. 



ratio of molecular to atomic gas (£h 2 /Ehi) is determined 
by the interstellar hydrostatic pressure. In order to com- 
pare Eh 2 /Shi with the hydrostatic pressure at z = 
(-Po)j we use the following approximation for Pq as de- 
rived in the Appendix [Bj 



to a factor of 5. The dashed lines in the figure will be 
discussed in Section HT4l 

A plot of Eh 2 /£hi against Pq (using constant o~ g and 
z*) shown in Figure [TTT b) suggests a power law relation- 
ship. We have fit a power-law slope of: 



P = 0.89(GE») U - 5 E 



S as 7 0.5 



(13) 



Ptr 



0.89 ± 0.01 



(14) 



where E* is the stellar mass density, o~ g is the gas velocity 
dispersion, and z* is the stellar scale height. This approx- 
imatio n closely matches that used bv lBlitz fe Rosolowskvl 
( 2004) to describe a two-component disk of gas and stars 
where the mass is dominated by stars having a thicker 
vertical distribution. 

We first assume in Equation [13] that <x g (= 8 km s _1 ) 
and z* (= 8") are c o nstan t with radius, as assumed by 
iBlitz fc Rosolowskvl (|2004[) . Observations suggest that 
stellar disks tend to have nearl y constant scale height 
rather than velocity dispersion: Ivan der Kruit fc Searld 
(|1981af ) found that the stellar scale heights in the edge- 
on spiral galaxies NGC 4244 and 5907 are largely inde- 
pendent of radius. We use the value of z*=8" that we 
obtained by fitting the exponential disk model in Section 
14.21 From the derived E gas and E* in Section [U we have 
obtained the radial pressure profile, normalized by the 
optical radius r25, shown in Figure [TlT a) as a red solid 
line (Model CG-ZS). We use the dash mark (-) for the 
combination of models. On the other h a nd, more recent 
studie s (e.g., Ide Griis fc Peletierl 119971 ; iNaravan fc Jod 
l2002a|) have suggested that the stellar scale heig ht in- 
creases with radius, and our analysis in Sect ion \h . 2 1 shows 
an increasing stellar scale height as a function of radius. 
Therefore, we also show the pressure when the gas ve- 
locity dispersion (<j g ) and the stellar scale height (z*) 
are varying with radius (obtained in Section [3]) in Figure 
rnja) as a black solid line (Model PG-PS). This model 
is associated with considerably higher pressures, by up 



3 K) is the pressure at the tran- 
= Ehi- Since we are explicitly 



Eh 2 
Ehi 

where P tl (~ 7xl0 3 cm" 
sition radius where Eh 2 
taking azimuthal averages, the uncertainties in R mo i are 
based on the standard deviation of the mean rather than 
the standard deviation of values in an annulus. The fit 
has been obtained with weights inversely proportional 
to the uncertainties shown in Figure [TTJb). The slope 
a = 0.89 of the weighte d fit is in reasonable agreement 
with previous studies bv IWong fc Bli tz (2002) \a = 0.8] 
IBlitz fc Rosolowskvl (|2006l ) [a = 0.92]^ and lLerov et al.l 
(2008) [a = 0.73]. In additio n, the value of P tI is within 
the range of values found bv lBlitz fc Rosolowskvl (2006) 
in their sample, although the value is about a factor of 5 
lower than the average value. 

6.2. Star Formation Rate and Efficiency 

In order to examine systematically the relationships 
between star formation rate and molecular, atomic, or 
total gas, we plot Esfr against Eh 2 , Shi and E gas at 
the same radius (using 12" wide radial bins) in Figure 
1121 showing that the correlation between Esfr and Eh 2 
is much stronger than that between EgpR and Ehi- The 
error bars in the upper left corner represent the uncer- 
tainties of the SFR and the gas surface densities, which 
are a factor of ~ 2. We have determined the Schmidt 
law index by least-squares fitting (solid line in Fig. [13 



to the data points of Esfr against S g as in order t o com - 
pare with the index (1.4) derived bv iKennicuttl (|1998|) . 
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Fig. 12. — SFR surface density as a function of H I (open circles), 
H2 (filled circles) and total gas (filled squares) surface density. The 
dashed lines represent constant SFE with the corresponding star 
formation timescale (1/SFE) labeled. The solid line shows the 
Schmidt law with index of 0.85: £ S FR oc E gas ' 85 ± °' 55 . The rms 
dispersion around the fit line is ~ 0.23 dex. The error bars in the 
upper left corner represent the uncertainties of the SFR and the 
gas surface densities. 



The Schmidt law index we have derived for the total gas 
is ~ 0.85 ± 0.55, while the index for the molecular gas 
is ~ 0.77 ± 0.44. The estimated scatter around the rela- 
tion is ~ 0.23 dex for the total gas and ~ 0.17 dex for the 
molecular gas. These sc a tter v alues agree well with those 
derived by iBigiel et al.l (|2008[ ). The gas depletion time 
for the total gas is in a range of 2.7-4.7 Gyr when both 
H2 and H I are included while the gas depletion time is 
in a range of 1.1-3.3 Gyr when only H2 is included. 

A simple theoretical estimate of the star formation rate 
is that it is proportional to the gas surface density T, g 
divided by the Jeans time, i.e., the time scale for the 
growth of gravitational instabilities. For a gas-only disk, 
the Jeans time is tj g ~ a g / '(7rC?E ff ), so 



^SFR,modl OC 



(15) 



where fb jff is the hydrostatic midplane pressure of the 
gas-only disk (see Appendix [Bj . As shown in Figure 
T37a). the correlation between the derived radial profile 
of Esfr and this theoretical estimate is relatively poor, 
regardless of whether a g is constant (Model CG) or vary- 
ing (Model PG), though a linear relationship with pro- 
portionality constant e=0.025 bisects the distribution of 
points. The poor correlation reflects the large scatter 
found in the observed Schmidt law. The representative 
error bars shown in Figure Q2] indicate a factor of 2 for 
Ssfr and a factor of 4 (from E|) for the theoretical es- 
timate. 

A more refined model considers the effe ct of both ga s 
and stars in calculating the Jeans time (jWongl (2009). 
This leads to a Jeans time of 



t.i. 



,r<7 [ ^ + 



(16) 



SsFR,mod2 OC 



c J,sg 



ttGIL 



Po 



(17) 



where Pq is the hydrostatic midplane pressure defined 
by Eq. (|B5[) . The correlation between the observed and 
predicted values of Esfr is now much better, as shown 
in Figure fT3Tb) . Indeed, the RMS difference (in the log) 
between the observed and predicted SFR has decreased 
from 0.100 dex to 0.005 dex. However, the predicted 
SFR is still much higher than the observed SFR in the 
inner disk (r < 120"), although the agreement is better 
when <jg is allowed to vary. The deviation from a linear 
relationship in t he inner dis k has been found in other 
galaxies as well ()Wong|[2~009f ). and reflects the fact that 
observed star formation time scales (1/SFE) are nearly 
constant in the inner, ^-dominated disks of galaxies, 
whereas dynamical timescales should become shorter as 
r decreases. We discuss this deviation further in Section 
O 

The SFE (SsFR/S gas ) as a function of radius is shown 
in Figure fTlT a). The SFE profile in the inner disk re- 
gion shows large variations, but these are due in part 
to uncertainties in the methods used to derive the £ gas 
and Ssfr profiles. The solid red profile shows the SFE 
using the total gas (£ gas ) profile from the PVD method 
and the dashed blue line presents the SFE obtained from 
using the RADPROF task to derive £ gas . They both 
show an overall decline in SFE in the outer disk. The 
vertical error bar in the upper right corner is obtained 
using the uncertainties of Ssfr and S gas , which are a 
factor of ~2. Despite differences up to a factor of ^4 
between the methods, the RMS difference is consistent 
with our adopted uncertainty. The dotted vertical line 
in the figure repr esents the tr a nsitio n radius (~ 240") 
where £h 2 = Shi- iLerov et alj (|2008l ) found, for a sam- 
ple of 12 spiral galaxies, an almost constant SFE inside 
the transition radius and decreasing SFE outside the ra- 
dius. Since SFR is linearly correlated with £h 2 (consis- 
tent with constant SFE for H2), the SFE is proportional 
to i?moi/(^?moi+l) and since i? mo i is a strong function of 
radius, SFE also depends on radius. The long-dashed line 
in the figure presents the best-fit relationship between 
SFE, R mn ] ( =Sho /Sht), and radius for spirals given by 
ILerov et all (|2~008): 



SFE = 5.25 x 10" 



-^mol 



R 



mol 



yr 



where 



Rmo\ = 10.6 exp(-r ga i/0.21r 25 ). 



(18) 



(19) 



Here r ga i is the galactocentric radius, r-^ is the opti- 
cal radius (~400" for NGC 891), and the expression for 
the SFE is based on the assumption of constant SFE in 
molecular gas. Within the errors, our data are consistent 
with the prediction. The predicted relationship between 
R mo \ and radius (Eq. HI?)) is shown in Figure H4T b) as 
the dot-dashed line alongside the Eh 2 /Shi ratio obtained 
from the PVD method (solid line) and the RADPROF 
task (dashed line). Our fitted exp onential scale lengt h 
of 80" is in good agreement with ILerov et"al] (2008), 
who found that i? mo i oc with a radial scale length 
of0.21r 25 . 
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Fig. 13. — (a) SFR surface density as a function of Po gaa — only/ a g (=irGS gas /2u 9 ) using Model CG (red squares) and Model PG (solid 
circles). The dotted line presents Ssfr proportional to Pogas-only/^g with efficiency 0.025. Note that the slope of the line is unity and 
the intercept indicates the efficiency. The representative error bars in the lower right c orner show a factor of 2 for S sfr and a factor of 
4 for the theoretical estimate, (b) Esfr as a function of the midplane pressure (Eq. IB5t divided by the gas velocity dispersion when 
Model CG and Model ZS arc adopted (red squares) and Model PG and Model PS are used (solid circles). The dotted line presents Ssfr 
proportional to Po/ctg with efficiency 0.012. 
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Fig. 14. — (a) Star formation efficiency (SFE) as a function of radius, SFE = SgpFi/Sg as . The solid red line shows the SFE as derived 
from S gas using the PVD method and the dashed blue line is the SFE as derived from S gas obtained from the RADPROF task. The 
dotted vertical line indicates the transition radius rt (~ 240"), where Ejj 2 equals Ejji in the PVD method. The dot-dashed line shows 
the predicted fit relationship between SFE and i? mo i given by Equation {l8j) . The representative error bar of the SFE is shown in the 
upper right corner, (b) Ratio of molecular to atomic gas surface density as a function of radius using the PVD method (solid line) and 
the RADPROF task (dashed line). The horizontal and vertical dotted lines show _R mo i = 1 and the transition radius, respectively. The 
dot-dashed line shows the predicted relationship between i? mo i and radius given by Eq. (1 1 9 D . 



6.3. Gravitational Instability 

Uog &: Solomonl |l984) have studied gravitational in- 
stability with two components consisting ol gas and stars 
(but both are collisional) in a galact ic disk t o deriv e an 
instability criterion. M ore recently. | Rafikovl (I2001D has 
extended the study of iJog &: Solomonl ( 19841) in order 
to investigate the instability criterion under axisymmet- 
ric gravitational perturbations in a more realistic galaxy, 
considering collisional gas and collisionless stars. We 



have derived the instability parame ter Q ga s-t -s tar us ing the 
following formulation provided by Kalikov (2001). The 
instability criterion for a thin rotating disk consisting of 
gas and stars is given by: 



-R(T 
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1 + q 2 Rl 



~[1 <T* 1 (q 2 )} > 1, (20) 
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Fig. 15. — (a) Q gas (blue dotted line with open circles), Qstar (red dashed line with asterisks) and Q gas +star (solid line with filled circles) 
radial profiles with a constant velocity dispersion of gas (Model CG) and stars (Model CS). The stellar velocity dispersion was calculated 
assuming cr* >r = <r„ jZ /0.6. Below the dotted line (Q = 1) lies the unstable region, (b) Q ga s, Qatar and Q gas +star radial profiles when 
Model CG and Model ZS are used, (c) Qgas, Qstar and Q gas +star radial profiles when Model PG and Model PS are employed. 



where 



9 



na 



7r6?S B 




k is the wavenumber (2tt/X), k is the epicyclic frequency, 
V is the rotational velocity, r is the galactocentric ra- 
dius, the radial stellar velocity dispersion, and Iq 
is the Bessel function of order 0. Note that the stellar 
velocity dispersion here is not in vertical direction but in 
radial direction. We assume this is relat ed to the vertica l 
velocity dispersion as a*. z = 0.6 x a* r (jBottemal 11991) . 
The rotation curve values shown in Figure [3] are used for 
the circular velocity V. We derive Q gas +star by choosing 
a specific A that yi elds a minimum value, following pre- 
vious studies (e.g., EI et al~ll200a lYang et al.ll2007[ ). All 
Q parameter profiles for gas, stars, and combination of 
gas and stars are shown in Figure [T5J We also compared 
Q profiles using a flat rotation curve (V = 250 km s _1 ) 
with the Q profiles shown in order to examine how much 
the inner region (within 100") Q profiles are affected by 
the rotation curve. The comparison shows that devia- 
tions from a flat rotation curve in the inner disk do not 
significantly affect the Q profiles. 

When constant velocity dispersions of gas (Model CG) 
and stars (Model CS) are adopted as explained in Section 
15.31 the Qgas+star radial profile is as presented in Figure 
[T5Ta). We also show the radial profile when Model CG- 
ZS and Model PG-PS are used in FigurefT5fb) and Figure 
fTBTc). respectively. The profiles in all the three cases 
show Qgas+star decreasing with radius in the inner disk 
and increasing in the outer disk. Note that regions where 



Q 



gas+star 



< 1 are unstable. We discuss our results in 



Section O 

6.4. Vertical Dependence of Gas Pressure 
The interstellar pressure as a function of (r, z), assum- 



ing turbulent support of the gas, is given by: 
P g (r,z) = p g (r,z)a 2 



where p g (r, z) is the gas volume density. The gas density 
is derived by summing the H2 and H I density profiles 
(pn 2 an d Phi, respectively), which are assumed to each 
follow a Gaussian distribution in z: 



p g (r, z) = 1.36 x 



(22) 



Po,h 2 M exp 



2/4 (r) 



Po.m{r) exp - 



2h 2 m (r) 



(21) 



Here the factor of 1.36 is a correction for helium, /5o,h 2 ( r ) 
and po,Hi('') are the midplane densities of H2 and H I 
function of radius, respectively, and h^ 2 (r) and 
hm (r) are the Gaussian widths of H2 and H I derived in 
Section 15.21 The densities at the midplane are derived 
from the radial gas surface density distributions obtained 
in Section|4]and from the disk thicknesses given in Figure 

Ha), using p = E/(W5tt). 

Figure [TH] shows the vertical profiles of P g at different 
radii when Model CG and Model PG are applied and the 
vertical profiles of the ph 2 /Phi ratio. The P g profiles 
show a pattern of values decreasing with radius near the 
midplane but increasing with radius at high z, regardless 
of model. Since P g oc p g for constant a g , this reflects the 
flaring of the gas disk at larger radii. R mo i (ph 2 /phi) 
in Figure 116( c) shows a similar pattern to the P g pro- 
files. However, the transitional z, where values change 
their pattern from decreasing to increasing with radius, 
is higher than that for P g , especially for Model CG. In ad- 
dition, values of i? mo i at high z show much less variation 
with radius than values of P g do. Green star symbols in 
Figure fTBT c) show where the vertically integrated values 
of i? mo i lie at each radius; the integrated value decreases 
with radius as seen in Figure H~4T b) . Figure [T7] shows 
-Rmoi( r j z ) against P g (r, z) at different heights using (a) 
Model CG and (b) Model PG. We focus on the region 
marked with filled symbols where the PVD and RAD- 
PROF methods give similar results for the radial profile. 
While the pressure using Model PG seems to govern R mo i 
over the range < \z\ < 10", the pressure using Model 
CG correlates with the ratio near the midplane but not 
necessarily at high z. 

The turbulent gas pressure at the midplane P g (r, 0)(= 
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Fig. 16. — (a) Dependence of gas pressure (using Model CG) on z at different radii from 90" to 210" in steps of 40". The horizontal 
dashed line represents the transitional z (~ 8"). (b) Dependence of gas pressure (using Model PG) on z at different radii from 90" to 210" 
in steps of 40". The transitional z is ~11". (c) Ratio of molecular to atomic volume gas density as a function of z at different radii from 
90" to 210" in steps of 40". The tr ansi tional z is about 16". The star symbols mark the ratio of molecular to atomic gas surface density 
at each radius obtained from Figure [I4jb) . 



pogCTg) based on Equation (f2T|) is shown as a function 
of radius in Figure lllf a) in comparison with the hydro- 
static pressure obtained in Section 16.11 using Equation 
(|13j) . When a constant value for o~ g is assumed, although 
the curves are similar in shape, there is a discrepancy in 
that Equation ([13]) predicts a factor of 2-3 larger value 
for Pq. On the other hand, the curves for radially varying 
a g (Model PG) are closer to each other, especially in the 
inner region. We discuss these trends further in Section 

7. DISCUSSION 
7.1. SFR and Midplane Pressure 

Our results in Figure Q~3] suggests that much of the 
scatter in the observed Schmidt law reflects real varia- 
tions in the relationship between Ssfr and £ g , and not 
just uncertainties in the determination of radial profiles 
for an edge-on galaxy. The time scale for star forma- 
tion depends not only on the properties of the gas disk 
but on those of the stell ar disk as w e ll, a point which 
has b een emphasized by iLerov et~aLl (|2008T ) and iWond 
(2009|). The basic reason that the midplane hydrostatic 
pressure correlates with the large-scale SFR is that both 
quantities reflect how quickly gravitational disturbances 
grow in a pressure-supported disk. 

How then can we explain the deviation from a linear 
relationship between SFR and Po/cr g seen in the inner 
disk (Figure H3b)? One possibility, based on the results 
shown in Figure [T5l is that an increase in the Q parame- 
ter towards the galaxy center suppresses gravitational in- 
stabilities, leading to a SFR that is still high but not quite 
as high as would be expected fro m the J e ans ti me scale. 
Another possibility, discussed by iWongj (|2009ft , is that 
the Jeans time scale is more relevant for the formation 
of giant molecular clouds (GMCs) rather than stars, and 
thus in the inner, H2-dominated disk the star formation 
time scale reflects other physical processes. We caution, 
however, that there are likely to be large uncertainties 
in £* and £ s in the inner disk, due to stellar population 
gradients and variations in the CO-to-H2 conversion fac- 
tor, and if these quantities have been overestimated by 



our simple linear conversion factors then the Jeans time 
scale will be underestimated, bringing the observed and 
predicted SFRs into better agreement. 

7.2. Q in the Star Formation Disk 

As shown in Figure [15j the Q gas +star radial curve ap- 
pears to increase toward the galaxy center, which means 
star formation should be suppressed in the center. There- 
fore, the Q curve shown in the figure seems not to match 
the SFR radial profile in Figure EJc), which shows active 
star formation in the center. This tendency, a increasing 
Q pro file toward the center, is also shown in ILerov et al.l 
(2008) who employ a constant stellar scale height and a 
stellar velocity dispersion depending on the square root 
of the stellar surface density (corresponding to Model 
ZS). This inconsistency between the <3 ga s+star and the 
SFR radial profiles is mainly due to the Q parameter 
depending on the epicyclic frequency k, which is propor- 
tional to r _1 for a flat rotation curve. Thus, k will tend 
to increase more quickly toward the center than the ra- 
tios of £/<t which govern the other components of Q. In 
addition, the presence of the stellar bulge could make the 
Q formulation, which assumes that matter is distributed 
in a disk, incomplete in the central region. 

The Qgas+star depends quite sensitively on the as- 
sumed velocity dispersions as seen in Figure 1151 The 
Q curves using Model CG-CS and Model CG-ZS show 
mostly values less than 1, meaning gravitationally un- 
stable, but their unstable regions are not the same. On 
the other hand, the Q values with Model PG-PS are 
mostly marginal (Q ~ 1) in the disk. Therefore, the 
Q profile using Model PG-PS where the velocity dis- 
persion is determined from the disk thickness appears 
cons istent with the idea of a self-regulating Q parame- 
ter (jKovama fc Ostrikeil I2009D . where high gas surface 
density and star formation in the central region leads to 
higher a g . 

7.3. Interstellar Gas Pressure in Two Dimensions 

The main difference between the turbulent gas pressure 
P g (r, 0) (=Po g o~g) and the hydrostatic midplane pressure 
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Fig. 17. — (a) Ratio of molecular to atomic volume gas density with the pressure (P 9 /k) using Model CG at different heights; red solid 
line at the midplane (z = 0"), blue dashed line at z = 5", and green dot-dashed line at z = 10". The triangle, square and circle symbols 
ind icate values at the radius 90", 150" and 210", respectively. The dotted line is the fitted line, proportional to p°- 89 , shown in Figure 
lllf b). (b) Ratio of molecular to atomic volume gas density with the pressure (P 3 /k) using Model PG at different heights. 



in Model CG shown in Figure flTTa) is that Equation 
(fl"3|) implicitly uses the stellar disk mass to determine 
the thickness of the gas disk, whereas with Equation (|23|) 
we determine the thickness of the gas disk directly. At 
R = 100", for example, the stellar surface density of 200 
M pc -2 and the adopted stellar scale height of 350 pc 
imply a stellar density of 0.3 M a pc -3 , which in turn im- 
plies a Gaussian width for the gas of 60 pc, or only ~1", 
which is a factor of 3—5 less than what we measure. To 
reconcile the two approaches, we must assume that we 
have (1) overestimated the stellar mass, by adopting too 
large a M/L ratio in Section 14.21 or (2) overestimated 
the gas thickness as a result of line-of-sight projection 
effects. (The other possibility, that we have underes- 
timated the stellar disk thickness, seems unlikely given 
that projection effects would tend to work in the oppo- 
site direction) . Future studies of edge-on galaxies should 
be aimed at testing these possibilities. 

Alternatively, the assumed constant a g in Model CG 
m ay be too low, as sugge sted also by the recent work 
by iTamburro" etafl (l2009h . We compare the pressures 
using a varying velocity dispersion o g (r) (Model PG) in 
Figure flTT a). There is still a discrepancy between the 
hydrostatic and turbulent gas pressures (black solid and 
dashed lines, respectively) in the outer region, but they 
are nearly equal in the inner region. The discrepancy is 
because the approximation for hydrostatic pressure (Eq. 
I13j) ignores the gas volume density term, which is im- 
portant in the outer disk. The profile shown as filled 
circles in Figure ITTT a') represents the hydrostatic pres- 
sure including the gas density (see Appendix Eq. IB6|) 
and is consistent with the turbulent gas pressure (black 
dashed line) , as it should be given that our derivation of 
<j g assumes hydrostatic equilibrium with only turbulent 
support. The fact that P g is very sensitive to the adopted 
<j g implies that comparison of ISM pressures deduced for 
face-on and edge-on galaxies must be made with caution. 
Fortunately, however, the predicted SFR in Eq. (fT"7| de- 
pends on P/cr g and not just P. 



The vertical gas pressure profiles at several radii shown 
in FigurefTBTa) and (b) are consistent with a recen t study 
focused on our Galaxy (Ka sparova fe; Zasovll2008l ) , which 
showed that the pressure decreases with increasing z at 
all radii, but that the sign of the radial pressure gradient 
reverses sign at high z (z > 8" for model CG and z > 
11" for Model PG; see left and middle panels of Figure 
[T6|) . due to flaring of the gas disk. In other words, the 
pressure decreases with radius when z is low, while the 
pressure increases with radius at high z. As shown in Fig- 
ure !17f a). which shows the correlation between i? mo i and 
P at various heights above the midplane, the pressure 
using Model CG (constant velocity dispersion) is corre- 
lated with i? mo i near the midplane but their relationship 
is not clear when z is large. However, the pressure using 
Model PG (varying velocity dispersion) in Figure flTT b) 
behaves more similarly to i? mo i- This is not surprising, 
since Model PG is based on the observed gas thickness 
used to derive the volume densities for i? mo i, although no 
distinction was made between CO and H I when deriving 
<J g {r). Since it is possible that UV radiation and metal- 
licity as well as the hydrostatic pressure affect the value 
of i? mo i, they could contribute to the relatively uniform 
vertical gradients in i? mo i seen at different radii. This 
may weaken the correlation between the gas pressure and 
R mo \ at high z. 

8. SUMMARY AND CONCLUSIONS 

We have derived the azimuthally averaged surface den- 
sity profiles for the CO, H I, and IR (3.6 and 24 /im) 
emission and the rotation curve to study the relationship 
between ISM and star formation in the edge-on galaxy 
NGC 891. In addition, we have estimated the gas volume 
density profile in two dimensions (r, z) using the mea- 
sured disk thickness, and inferred velocity dispersions as 
a function of radius in this galaxy. 

1. We have explored the vertical structure by fitting 
single or double Gaussian profiles to the CO and H I 
maps. The integrated H I data have been fitted by 
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a double Gaussian profile, implying two components of 
thin and thick disks. On the other hand, the CO disk 
has only one component: a thin disk, although sensitive 
single-dish mapping is still needed to confirm this. 

2. We have investigated the relationship between the 
interstellar hydrostatic pressure in the midplane and the 
ratio of molecular to atomic surface mass density and 
found a power law relationship with slope a — 0.89. 

3. The SFR surface density and molecular surface 
density profiles show similar behavior. In addition, the 
plot showing relationships between Esfr and £h 2 > ^hi, 
and Sgas presents a strong correlation between Ssfr and 
Sh 2 - The Schmidt law index we obtained using E gas is 
0.85±0.55, but there is considerable scatter around this 
relation. The power l aw index may be smaller than the 
index (1.4) derived by lKennicuttl (|1998f ) because we mea- 
sure the slope over mostly H2 dominated regions. 

4. The SFR surface density is more strongly correlated 
with the hydrostatic midplane pressure including both 
gas and stars and using varying velocity dispersions of gas 
(Model PG) and stars (Model PS), but the correlation 
appears to break down in the inner region. 

5. In order to study how the instability parameter 
(Qgas+star is related with star formation disk, we have 
derived Q ga s+star in a thin rotating disk model consist- 
ing of both gas and stars. The Q ga s+star radial profile 



with varying velocity dispersions for the gas and stars 
is favored from the point of view of leading to marginal 
instability throughout the disk, but still predicts a sup- 
pression of star formation near the center which is not 
apparent. 

6. The ph 2 /phi ratio against turbulent gas pressure 
using constant velocity dispersion (Model CG) has been 
compared with that using varying velocity dispersion 
(Model PG) in Figure [i~7l The latter model with varying 
o g correlates much better with the ph 2 / Phi ratio over 
the range < \z\ < 10". 

7. Estimates of the hydrostatic midplane pressure Pq 
based on a constant a g appear to substantially underes- 
timate the actual turbulent pressure needed to explain 
the thickness of the gas disk. This implies that either 
tjg is significantly higher than usually assumed, or that 
other sources of support for the gas disk are important. 
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APPENDIX 

RADIAL DISTRIBUTIONS OF H 2 AND H I 

Comparison between PVD and RADPROF 

In addition to the radial gas profiles derived by the PVD method shown in Fig. Eta), we obtained the radial 
distribution of H2 and H I using the RADPROF program to compare the two methods. Figure [18] shows the comparison 
between the two methods (PVD and RADPROF). The RADPROF profiles appear roughly similar to the radial profiles 
obtained from the PVD method. But there is a discrepancy in the central region for H I, where both methods are 
less reliable. RADPROF suffe rs from smoothing effects and is unable to accurately reproduce small-scale structure, 
as discussed in Warmcls] (|1988l ). Also, it works from the outside inward and errors can accumulate in the center. Note 



that S/N of H I is low at the center in contrast to the H2 profile. The PVD method is also less reliable in the central 
region, since it relies on the assumption of a flat rotation curve, which is likely to be incorrect in this region. Also, since 
a fixed velocity width corresponds to a smaller line-of-sight depth near the center of the galaxy, especially near the 
terminal velocity, the face-on brightness of a PVD pixel becomes higher if located near the center of the galaxy. This 
will tend to magnify noise fluctuations near the center. However, total flux of H I appears to be not much different 
(within ~ 5%) for the two methods in spite of the discrepancy near the center. We estimate that the uncertainty in 
the total gas profile is a factor of 2 based on the standard deviation of the differences between PVD and RADPROF 
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Fig. 18. — (a) Comparison between the RADPROF and PVD methods for H2 radial profile, (b) Comparison between the RADPROF 
and PVD methods for H I radial profile. 
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profiles of CO shown in Fig. ITST a). Since H I does not contribute much to the total gas, we have estimated the 
uncertainty from only the CO profile. We classify this uncertainty as a potential systematic error in the PVD method, 
to be distinguished from errors due to thermal noise or non-axisymmetry of the galaxy. 

Galaxy Models 

We have generated models of NGC 891 using the GIPSY task GALMOD to examine how well the two methods 
for deriving radial profiles are able to recover an input model. CO and H I data cubes are used as an input to set a 
coordinate system for the CO and H I models. Each model consists of rings located at every 10" from the center up 
to 280" and each ring is formed by a circular velocity, a surface density, and several properties such as scale height, 
inclination, velocity dispersion, and position angle. For the surface density we used the radial profiles obtained by the 
PVD method in Figure [UJa). The circular velocity at each ring is the assumed value of 250 km s" 1 , which is used in 
the PVD method. For the other properties such as scale height (19" for HI and 7" for CO), inclination (90°), velocity 
dispersion (8 km s^ 1 ), and position angle (23°), we used a constant value at each ring. Figure fT9l shows radial profiles 
of the models (derived from the two methods PVD and RADPROF) and the input surface-density profile. In general, 
the retrieved PVD profiles from the CO and H I models match their input profiles. Also, the profiles obtained by the 
two different methods agree each other well, although the RADPROF profile is more smooth because RADPROF does 
not use velocity information. 
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Fig. 19. — (a) Radial profiles of CO model obtained from the PVD method (blue line with open circles) and the RADPROF task (red 
line with open squares). The dashed line profile with solid circles represents the input profile for the CO model, (b) Radial profiles of H I 
model obtained from the PVD method (blue line with open circles) and the RADPROF task (red line with open squares). The dashed line 
profile with solid circles represents the input profile for the H I model. 



APPROXIMATIONS FOR HYDROSTATIC PRESSURE 



Close to th e midplane of a self-gravitating disk with an isothermal velocity dispersion a, the density is given by 
(|Spitzerlll942l) 



p = Po sech — 

\*0 



where 



u 2nGp 

Integrating in z yields a mass surface density E = 2pqZq. The turbulent pressure at the midplane is then given by 



(Bl) 



c 2 2 

Po = Pqo- = — E 



(B2) 



The case of a gas disk embedded in a stellar disk (jTalbot &: Arnettlll975| ) modifies the density distribution to be 

>(z)-<X0) 



Pg( z ) = Pog exp 



(B3) 



where <fi is the gravitational potential. Taking the leading term in Poisson's equation when expanding around z = 
yields a Gaussian distribution 
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Note that in this case the g as surface density is S g 
gas pressure to be (see also Elmcgrcenl ll989ft : 



(2TT)°- 5 p ghg. iTalbot fc Arnetti (|1975|) then derive the midplane 



ttG 



ttG. 



(B5) 



where the expression in square brackets is approximately the vertical oscillation frequency of a test particle in the 
combined potential of gas and stars. If the dominant stellar disk maintains a sech 2 density distribution (Equation [Bl}, 
this expression can be recast in terms of densities as 



Pn = 



ttG 



Pog 



2(p *) 



0.5 



L(2Gp ,- 



\0.5 



(2ttG) -' 



(B6) 



In the star-forming disk of a large spiral galaxy, neglecting the influence of dark matter in the disk, one can often 
assume po g -C po* ~ Po.tot- Thus, ignoring the first term in brackets leads to 



%G 



^g a g 



(B7) 



which is approximately the expression adopted by iBlitz fc Rosolowskvl ()2004l ) . Note that iBlitz fc Rosolowskvl (|2004fi 
define their stellar scale height as h r = z^j\[2. Note also that a g should technically be co nsidered the "effective" 
velocity dispersi on to allow for partial support from cosmic-ray and magnetic field pr essure ([Talbot fc Arnetti 119751 : 
lElmeg rccn 1989), which can increase the gas scale height. However, [Elmcgrccn (|1989l) has argued that the external 
pressure on an interstellar cloud is primarily the kinematic component, which is of order 60% of Pq. 
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